ON THE SOLUTIONS OF WEAK NORMALITY 
EQUATIONS IN MULTIDIMENSIONAL CASE. 



RUSLAN A. SHARIPOV 

Abstract. The system of weak normality equations constitutes a part in the com- 
plete system of normality equations. Solutions of each of these two systems of equa- 
tions are associated with some definite classes of Newtonian dynamical systems in 
Riemannian manifolds. In this paper for the case of simplest flat Riemannian mani- 
fold M = R n with n ^ 3 we show that there exist solutions of weak normality equa- 
tions that do not solve complete system of normality equations in whole. Hence 
associated classes of Newtonian dynamical systems do not coincide with each other. 



1. Introduction. 

Let M be a Riemannian manifold of the dimension n, and let S be a hypersurface 
in M. One of the ways for deforming S consists in moving points of S along 
trajectories of some Newtonian dynamical system. In local coordinates x , ... , x n 
in M such system is given by n ordinary differential equations 

n n 

(1.1) ■ XX 'V^' F V >■"■■ 

i=i j=i 

where k = 1, . . . , n. Here = r^x 1 , . . . , x n ) are components of metric connec- 
tion for basic metric g of the manifold M, and F k are components of force vector 
F. They determine force field of dynamical system (1.1). 

In order to obtain a shift of hypersurface S we set up the following Cauchy 
problem for the system of ordinary differential equations (1.1): 

(1.2) x k \ = x k {p), x k \ = v(p)-n k (p). 

lt=o lt=0 

Here p is a point of hypersurface S, while n (p) are components of unitary normal 
vector to S. The function v(p) is interpreted as modulus of initial velocity v|; we 
assume it to be a smooth function of p. For the fixed point p initial data (1.2) 
determine some trajectory of dynamical system (1.1) coming out from the point p 
and being perpendicular to S at this point. Let's map the point p to the point p(t) 
on such trajectory. If such correspondence can be extended to the whole surface S, 
then we have a shifting map ft : S — > St- In general case the theorem on smooth 
dependence of the solution of ODE upon initial data (see [1] or [2]) warranties 
only the possibility to extend this map to some neighborhood of the point p in 
S. Let S' be such neighborhood. When t is sufficiently close to zero, shifting 
maps ft: S —y St are diffeomorphisms, their images are smooth hypersurfaces. In 
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whole, diffeomorphisms f t : S' — > S' t , which are defined locally, are glued into a 
one-parametric family of local 1 diffeomorphisms ft'. S —> St- 

Definition 1.1. One-parametric set of local diffeomorphisms f t :S—>St deter- 
mined by the equations (1.1) and by initial data (1.2) is called a shift of S along 
trajectories of dynamical system (1.1). Such shift is called a normal shift if all 
hypersurfaces St produced by this shift are perpendicular to its trajectories. 

Let po be some point of hypersurface S, and let v$ be some nonzero number. We 
normalize the function v(p) in (1.2) by the condition: 



(1.3) 



Definition 1.2. Newtonian dynamical system (1.1) with force field F is called a 
system admitting the normal shift in strong sense 2 if for any hypersurface 
S in M, for any point po G S, and for any real number iv ^ there exists a 
neighborhood S' of the point S, and there exists smooth nonzero function v(p) in 
S' normalized by the condition (1.3) and such that the shift f t : S' — > S' t determined 
by this function is a normal shift in the sense of definition 1.1. 

Definitions 1.1 and 1.2 underlie in the base of the theory of dynamical systems 
admitting the normal shift. This theory was developed in papers [3-18]; results 
of these papers were used in preparing theses [19] and [20]. In paper [8] we have 
derived the following equations, which were called weak normality equations: 



i=l 



j=i 



(N*F j )]Pi=0, 



(1-4) { J2Y,(V i F j +V j F i -2v- 2 F i F j )WPZ + 

i=l j=l 

i=l j = l \ r=l / 

In paper [9] additional normality equations were derived: 

, n n / n ~ \ 



(1.5) 



t — 1 j — 1 \ 771—1 

n n / 



7 — 1 j—1 \ m= 



EE^v^pf = EEE 



F- V F 



v i=l j = l 



2—1 j — 1 m—1 



n - 1 



PI. 



1 Locality here means that domain of the map ft depends on t. For sufficiently large ( it can 
be empty at all. 

2 Earlier we used the definition without normalizing condition (1.3) for the function v(p). Such 
definition was called the normality condition. The definition 1.2 strengthens this condition making 
it more strict with respect to the choice of force field F of the dynamical system (1.1). It is called 
the strong normality condition. 
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Leaving the equations (1.4) and (1.5) with no comments for a while, now we shall 
formulate two theorems binding these equations with definitions 1.1 and 1.2. 

Theorem 1.1. Newtonian dynamical system in two-dimensional Riemannian ma- 
nifold M satisfies strong normality condition if and only if its force field F satisfies 
weak normality equations (1.4) for v = |v| ^ 0. 

Theorem 1.2. Newtonian dynamical system in Riemannian manifold M of the 
dimension n 3 satisfies strong normality condition if and only if its force field F 
satisfies the normality equations (1.4) and (1.5) for v = |v| ^ 0. 

Theorems 1.1 and 1.2 show that two and multidimensional cases are substantially 
different. In multidimensional case complete system of the equations (1.4) and (1.5) 
is strongly overdetermined, so that in reducing it we find that it is integrable in 
explicit form. General solution for the equations (1.4) and (1.5) is determined by 
two arbitrary functions h = h(w) and W — W(x 1 , ... ,x n ,v): 



Formula (1.6) was obtained in thesis [19]. Here ViW is the partial derivative of the 
function W in i-th coordinate of the point p, while W v is partial derivative of W 
in the variable v, which is is interpreted as modulus of velocity vector: v = |v|. 

In two-dimensional case we have only the equations (1.4). Here they were reduced 
to one scalar partial differential equation of the second order, then various special 
solutions of this equation were constructed (see thesis [20]). Most of these solutions 
correspond to the force fields that can be obtained by means of formula (1.6) when 
taking n = 2 in it. These force fields were called fields of multidimensional type. 
However, one of the most important results of thesis [20] is that there essentially 
two-dimensional solution of the equations (1.4) was constructed. This solution is 
not given by formula (1.6). Here we also construct the solution of the equations 
(1.4), which is not expressed by formula (1.6), but in multidimensional case. 

Formula (1.6) obtained in thesis [19] describes all Newtonian dynamical systems 
admitting the normal shift in Riemannian manifolds of the dimension n ^ 3. How- 
ever, when thesis [19] has been already written, new statement of the problem of 
normal shift was found, ft leads to the equations (1.4) in pure form (without ad- 
ditional equations (1.5)). Indeed, in papers [3-18] and in theses [19] and [20] we 
considered only smooth hypersurfaces S and we choosed sufficiently small values of 
t for the sift ft : S — > St to result in smooth hypersurfaces St only. If we eliminate 
this restriction for t, then in the process of shifting we sometimes can observe sin- 
gular points on hypersurface St (they arc called caustics). In particular, under the 
definite sircumstances hypersurface St can contract into a point at a time for some 
t = to- This process is called collapse. Immediately after the callapse for t > t 
we shall observe a blow-up of the point into a series of expanding hypersurfaces 
S t . The idea to consider the blow-ups of points by means of Newtonian dynam- 
ical systems were suggested to me by A. V. Bolsinov and A. T. Fomenko when I 
was reporting results of thesis [19] in the seminar at Moscow State University in 
February of 2000. This idea was realized in paper [21] and in paper [22]. 

Let po be some point of Riemannian manifold M. Let's consider the set of all 
unitary vectors in tangent space T po {M). They can be interpreted as radius- vectors 



(1.6) 




4 



R. A. SHARIPOV 



of the points of unit sphere a in T po (M). Let q £ a and let n(g) be the radius- vector 
of the point q on unit sphere. Let's fix some constant number v ^ and set up the 
following Cauchy problem for the equations of Newtonian dynamical system (1.1): 

(1.7) x k \ = x k {p a ), x k \ =v Q -n k (q) 

lt=0 lt=0 

For the fixed q initial data (1.2) determine some trajectory of dynamical system 
(1.1) coming out fromn the point pq. To the point q £ a we put into correspondence 
the point p(t) on such trajectory. Theorem on smooth dependence of the solution of 

ODE upon initial data (see [1] or [2]) says that 
we can extend this map to some neighborhood 
of the point q in a. If a' is such neighborhood 
of the point q, then we have one-parametric 
family of diffcomorphisms /( : u' — ► S' t . Due 
to compactness of unit sphere a we can glue 
local maps into one map f t : a — > S t . Here 
diffcomorphisms f t : a — * St are determined 
globally on the whole sphere a, though pa- 
rameter t can be restricted by some interval 
(— £, +e) on real axis as before. 

Definition 1.3. One parametric family of dif- 
feomorphisms f t : a — > S t given by the equa- 
tions (1.1) and initial data (1.7) is called a blow-up of the point p along trajecto- 
ries of dynamical system (1.1). It is called a normal blow-up if all hypersurfaces 
St arising in this blow-up are perpendicular to its trajectories. 

Definition 1.4. Newtonian dynamical system (1.1) with force field F in Ricmann- 
ian manifold M is called admitting normal blow-up of points if for any point 
Po £ M, and for any positive constant vq initial data (1.7) determine normal blow- 
up of this point along trajectories of dynamical system (1.1). 

Definitions 1.3 and 1.4 were first formulated in paper [21]. They introduced new 
object: a class of newtonian dynamical systems admitting the normal blow-up of 
points in Riemannian manifolds. In paper [21] it was shown (see theorem 12.1) 
that this new class of systems comprises the class of dynamical systems admitting 
the normal shift of hypersurfaces, which was previously considered. More exact 
description of new class of dynamical systems is given by the following theorem 
proved in paper [22]. 

Theorem 1.3. Newtonian dynamical system (1.1) on Riemannian manifold M 
admits normal blow-up of points if and only if its force field F satisfies weak nor- 
mality equations (1.4) for |v| ^ 0. 

From theorems 1.1 and 1.3 we see that class of dynamical systems admitting the 
normal shift of hypersurfaces in two dimensional case n = 2 coincides with the class 
of systems admitting normal blow-up of points. This was proved in [22]. But for 
the multidimensional case n ^ 3 the question on coinciding or not coinciding 
of these two classes remained open. I. A. Taimanov was strongly interested in this 
question during my report in the seminar of Yu. G. Reshetnyak at MI SB RAS 
(Mathematical Institute of Siberian Brunch of Russian Academy of Sciences) in 
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October of 2000. The main goal of present paper is to give answer to this question 
and, thus, eliminate one more obstacle for defending 1 thesis (19). 



Let's consider the force field F in the equations of Newtonian dynamics (1.1). 
Left hand side of these equations are components of acceleration vector VjV (covari- 
ant derivative of velocity vector with respect to parameter r along the trajectory). 
Therefore F k are components of tangent vector to M. However, they depend on 
double set of arguments: on coordinates x 1 , ... , x n of the point p S M and on 
the components of tangent vector v G T p {M). Pair q = (p, v) is a point of tangent 
bundle TM, so that p = n(q). Thus, considering the equations of the form (1.1), 
we come to the concept of extended vector field. Its generalization is a concept 
of extended tensor field. 

Definition 2.1. Function X that maps each point q = (p, v) of tangent bundle 
TM to some tensor of the type (r, s) from the space TJ(p, M) at the point p = ir(q) 
is called extended tensor field of the type (r, s) on the manifold M. 

The concept of extended tensor field stems from thesis [23] of Finsler, which 
gave rise to Finslerian geometry. Another approach to constructing extended ten- 
sor fields consists in rising them from M to TM. Here they constitute some special 
subset in the set of traditional tensor fields on TM. Such tensor fields were con- 
sidered in the book [24], they were called semibasic tensor fields. I am grateful 
to N. S. Dairbekov from IM SB RAS, who noted that theories of extended and 
semibasic tensor fields are isomorphic to each other. 

Below we shall use theory of extended tensor fields, which is described in details 
in Chapters II-IV of thesis [19]. It is based on the definition 2.1. 

For the beginning let's consider some particular examples of extended tensor 
fields on the Riemannian manifold M. 

1. Let's take the point q = (p, v) of tangent bundle TM and let's map it to the 
vector v belonging to tangent space T p (M). This yields an extended vector field, 
which is called the field of velocity. 

2. Let's map the point q = (p, v) of TM to the number v — |v|. This yields an 
extended scalar field, which is called the field of modulus of velocity vector. 

3. Extended field of unitary vectors N is determined as the ratio of two previous 
fields: N = v/w. 

4. Extended field of operators P is formed by operators of orthogonal projection 
onto the hyperplanes perpendicular to velocity vector v. Its components can be 
written explicitly: P) = 5) - N i Nj. 

Components of all above fields are present in normality equations (1.4) and (1.5). 
Moreover, in these equations we see the operators of covariant differentiation V and 
V; they called spatial and velocity gradients. The simplest way to define them is 
to use explicit formulas in coordinates: 



2. Normality equations and extended tensor fields. 



(2.1) 



VmX n'-X - 




dv 



1 The matter is that in Russia one can pretend for the degree of Doctor of Sciences only upon 
writing thesis and passing so called Defense Procedure in Specialized Council. 
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p,^ri\...% T n n flY 11 '- 

(2.2) Q=lb=1 

V / r n s n 

. pifc t^u... a fc ... i r pbfc v 

^ 2-^1 2-^1 ma k j s mjl 



fc=l «i=l fc=16fc=l 



3. Scalar ansatz. 

In order to simplify the normality equations (1.4) and (1.5) in paper [18] the 
scalar ansatz was suggested. Note that it is applicable either in two-dimensional 
case n = 2, and in multidimensional case n > 3 as well: 

n 

(3.1) F k =AN k -\v\^2PiViA. 

i=l 

Scalar ansatz (3.1) follows from first normality equation in the system (1.4). Sub- 
stituting (3.1) back into this equation, we turn it into identity. 

Formula (3.1) expresses force vector F through one extended scalar field A. This 
scalar field A can be expressed back through F: 

n 

(3.2) A = ^F i N i . 

i=l 

Substituting (3.1) into the second equation in the system (1.4), we obtain the 
following equation for scalar field A: 

n / n n 

( V.A + I v| ]T ]T P qr V q A \7 r \7 s A - 

(3.3) 5=1 V q=lr=1 

n n \ 

-^2N r A V r V s A - |v| Nr V r V S A P fe s = 0. 

r— 1 r— 1 / 

Formulas (3.1) and (3.2) establishes one-to-one correspondence between solutions of 
the equations (1.4) and (3.3). If force field F satisfies the equations (1.4) and (1.5) 
simultaneously then it is expressed by formula (1.6). Scalar field A corresponding 
to such force field is given by formula 

(3 - 4) A= ^r- V ^^— 

i=i 

Formula (3.4) admits gauge transformations, which change functions h = h(w) and 
W = W(x^~ , ... , x n , v), but which don't change A: 

W(x\... ,x n ,v) -^p(W(x\... ,x n ,v)), 

(3.5) 

h(w) — > h{p- l {w)) P \p- l (w)). 



If the function h(w) is nonzero, then by means of gauge transformations (3.5) it 
can be made identically equal to unity (see thesis [19] and succeeding paper [21]). 
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Therefore, instead of formula (3.4) with two arbitrary functions, we can use two 
formulas with one arbitrary function: 



(3.6) A = < 



1 " N i VtW 

i=l 



-v > — — for h = 0. 



Now the problem stated in section 1 is reformulated as follows: one should find a 
solution of the equations (3.3) that cannot be expressed by formula (3.6) neither 
for h = 1 nor h = 0. 

4. Spatially homogeneous force field with axial symmetry. 

Let M be a space R n with standard Euclidean metric. We shall construct the 
required solution of the equations (3.3) in this simples case. In the space R" co- 
variant derivatives V m and V m given by formulas (2.1) and (2.2) turn to partial 
derivatives: V m = d/dx m and V TO = d/dv m . We restrict our consideration to 
spatially homogeneous force fields, in R n they depend only on velocity vector, but 
don't depend on coordinates x 1 , ... , x n . For the corresponding function A this 
yields A = A(t> 1 , . . . , w"). When substituting this function into the equation (3.3), 
this equation is reduced to the following one: 

n / n n n \ 

(4- 1 ) E U V I E E ^ A VrV.A -J2 N " A VrV s A P S k - 0. 

s — 1 \ q— 1 r— 1 r— 1 / 

Let's mark some direction inIR™ determined by some constant unitary vector m. 
Without loss of generality we can assume m to be directed along n-th coordinate 

axis. Let's expand v into a sum 

(4.2) v = u ■ n + w ■ m, 

where n ± m and |n| = 1 All directions per- 
pendicular to m are assumed to be equiva- 
lent. Therefore we choose function A depend- 
ing only on two variables: A = A(u,w). One 
should study whether such choice is compat- 
ible with normality equations (4.1). For this 
purpose we calculate partial derivatives 



(4.3) V t A= dv 



dA \ A u ■ — for i < n 



u 

for i = n. 



Formula (4.3) shows that vector of velocity gradient VA belong to the linear span 
of vectors m and n from the expansion (4.2): 

(4.4) VA = A u ■ n + A w ■ m. 
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Denote by B the vector with the following components: 

n 
9=1 

This is the projection of velocity gradient \7A to the hyperplane perpendicular to 
velocity vector. One can obtain explicit formula for the vector B. For this purpose 
let's denote by 9 the angle between vectors m and v. Then let's expand projections 
Pm and Pn in the base composed by unitary vectors m and n: 

Pm = sin 2 9 ■ m — sin 9 cos 9 ■ n, 

(4.5) 

Pn = cos 9 ■ n — sin 9 cos 9 • m. 
From (4.4) and (4.5) we derive the following expression for B: 
B = (A u cos 2 9 — A w sin 9 cos 9) • n + 

(4.6) 

+ (A w sin 9 — A u sin# cos#) • m. 

Components of the vector B are present in the equation (4.1). Now we can write 
this equation in the following form: 

n n 

(4.7) Y^H^^ - AN ^ WF t s = 0. 

s—l r—1 

Now let's calculate the derivatives V r V s A, which are present in the equation (4.7). 
In order to do it we rewrite formula (4.3) as follows: 

V S A — A u h A w ui s . 

u 

Here (m |v) is the scalar product of vectors m and v. Now for second order deriva- 
tives V r V s ^4 by direct calculations we get 

~ v r - m r (m |v) v s - m s (m |v) 

U U 

, , v s -m s (m\v) v r — m r (m|v) 

A uw I 7Yl r -|- 77Z S I -|- A ww Ulr Tfl s . 



The above expression is rather complicated. In order to simplify it let's note that 
(m |v) = v cos 9, (n |v) = v ain9 — u. Hence 

v r - m r (m |v) _ n r (n |v) _ 



When applied to second order derivatives V r V s A, this equality yields: 
V r V s A = A uu n r n s -\- A uw {vn T n s -\- n r tti s ) -\- A ww in r m s . 
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Denote b = cos 9 • n — sin 9 ■ m. This is unitary vector belonging to the linear span 
of vectors v and m and being perpendicular to v. Using this vector we can rewrite 
the relationships (4.5) as follows: 

Pn = cos 9 ■ b, Pm = — sin 9 ■ b. 

Now we are able to contract \7 r \7 s A with the components of projector P: 

n 

V r V S A P k = A uu cos 9n r b k — A uw sin 9n r b k + 

(4.8) fi 

+ A uw cos 9m r b k — A ww sin 9m r b k . 
Let's use (4.8) for to rewrite (4.7) in more explicit form: 



( A uu cos 9 - A uw sin 9) (v B r - AN r ) n r b k + 

(4-9) 

+ y^( A uw cos 9 - A ww sin 9) (v B r - A N r ) m r b k = 0. 

r=l 

Vector b is nonzero. Its components cannot vanish simultaneously. Hence in left 
hand side of (4.9) we can collect common multiple b k and cancel it. Sums in r are 
scalar products. Therefore we get 

( A uu cos 9 - A uw sin 9) (v B - A N | n) + 
' " ' + (A uw cos6> - A ww sm9) (uB - AN | m) = 0. 

In order to calculate scalar products in (4.10) we use the expansion (4.6) for the 
vector B and the expansion N = cos 9 ■ m + sin# • n for the vector N: 

(A uu cos 9 — A uw sin 9) {v A u cos 2 9 — 

(4.11) —vA w sin6> cos# — A sa\9) + (A uw cos9 — A ww sin^)x 

x (v A w sin 2 9 — v A u sin 9 cos 9 — A cos 9) = 0. 

Theorem 4.1. System weak normality equations (3.3) written with respect to the 
function A = A(x , . . . , x n , v , . . . , v n ) in flat Euclidean case M — R" admits the 
substitution A = A(u, w), where u = ^/ (v 1 ) 2 + . . . + (v n ^ 1 ) 2 and w = v n . Thereby 
it is reduced to the single differential equation (4.11) ; where v — \J (u) 2 + (w) 2 and 
9 = arccos(w/w). 

The equation (4.11) is still rather complicated. In order to simplify it we trans- 
form it to polar coordinates v and 9 in the plane of variables u and w, i. e. we do 
the following change of variables 

(4.12) u = v sin#, w = v cos9. 

Let's calculate whether how partial derivatives are transformed under the change 
of variables (4.12). For the first order derivatives we have 

a a ' /\ a cos# sin 9 
A u = A v sm 9 + Ag , A w = A v cos 9 — Ag . 
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Then let's calculate second order partial derivatives: 

2 sin 20 cos 2 cos 2 sin 20 

Am = A™ Sin + h Agg - h Ay Ag , 

V V z V V 

. . sin 20 cos 20 sin 26* sin 20 cos 20 

Alto = A vv h A v g Agg — — A v — Ag — , 

2 v 2v z 2v v z 

a a in a sm $ a sin 2 sin 2 sin 20 

4» = A™ cos - Ae h Ae — 5— + A v h A — — . 

V v z v v z 

And finally, let's substitute all the above expressions for partial derivatives into the 
equation (4.11). Thereby the equation (4.11) crucially simplifies and takes the form 

(4.13) ^ + ^L+A e Ay = AAye. 

V V 

Theorem 4.2. System weak normality equations (3.3) written with respect to the 
function A = A(x x , . . . , x n , v 1 ,. .. , v n ) in flat Euclidean case M = K" admits the 
substitution A = A(v,6), where v = |v| and = arccos(u"/|v|). Thereby it is 
reduced to the single differential equation (4.13). 

Note that the equation (4.13) do not depend on the dimension of the space 
M = M™. It holds either in two-dimensional case n = 2, and in multidimensional 
case n ^ 3 as well. Moreover, this equation is well known in the theory of dynamical 
systems admitting the normal shift (see paper [10] and thesis [20]). In paper [10] 
was shown that the equation (4.13) is integrable in quadratures. For this purpose 
it was first transformed to the following form: 

Ag (A e \ _ v fAe\ + Ae + (Ae\ 3 = Q . 



A \Aj e \AJ V A \A 

This form of the equation (4.13) says that we should denote Ag/A = b. Then for 
b = b(v, 0) we obtain quasilinear partial differential equation of the first order: 

(4.14) bbg -vb v + b + b 3 = 0. 

In the equation (4.14) it is convenient to do another one change of variables, taking 
b = cotanz, where z — z(v, 0). This brings the equation (4.14) to the form 

(4.15) z e -v — z v -l = 0. 

cosz 

The equation (4.15) with the use of method of characteristics (see [25]). Charac- 
teristics of the equation (4.15) are the solutions of the following system of ODEs: 



(4.16) 



= 1, 

z = l, 



smz 
v = — V . 

cos z 
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Note, that from (4.16) one can derive the following equalities: 

■ . „ cos z . sinz . 
6-z = 0, —*-v-\ z = 0. 

IT V 

These equalities can be integrated. They mean that the system of equations (4.16) 
has the pair of first integrals I\ and I2 ■ 

(4.17) h=9-z, J 2 = — . 

v 

General solution of the equation (4.15) is determined by first integrals (4.17) in 
implicit form by means of functional equation 

(4.18) *(Ji,/ 2 )=0, 

where $ is some arbitrary function of two variables. 

Formula (4.18) proves the integrability of the equation (4.13) in quadratures. 
But in general it yields only local solution of this equation. Our goal is to construct 
global solution of the equation (4.13). It should be a smooth function in direct 
product of two intervals: closed interval [0, tt] for the variable 9 and open interval 
(0, +00) for the variable v. At the ends of the interval [0, n] one should provide 
the boundary conditions 

(4.19) A e \ =0, A e \ = 0. 

They appear because the function A(v, 9) should correspond to to the function of 
several variables A^ 1 ,... ,v n ) with axial symmetry being its restriction to the 
plane passing through the axis of symmetry. 

Now let's proceed with constructing the required solution of the equation (4.13). 
Note that functional relation of two first integrals I\ and I2 written as (4.18) for 
some cases can be given by the function of one variable y = f(w). Let's write (4.18) 
as I\ — fih), i- e. let's consider the following functional equation: 

,4.20) ,(22i). 

For f(w) we choose smooth increasing function with decreasing derivative; we as- 
sume that f(w) is defined in semiopen interval [0, +00) and /(0) = ir/2. Let's also 
assume that f(w) is not restricted and grows to infinity asun +00. Its graph 
is shown on Fig. 4.2. By means of function y = f(w) we construct a family of 
functions y = F [v] (z) depending on v as parameter: 

(4.21) y = F lv] (z) = f{^). 

Graphs of the functions are shown on Fig. 4.3. We use them in order to solve the 
equation (4.20) graphically. For this purpose we consider a family of straight lines 
being graphs of the following functions: 



(4.22) 



y = F m (z) = 9- z. 
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Denote by u min the value of the derivative f'(w) at the point w = 0: 

v^ = f(w)\ =0. 

\w=0 

Suppose that the values of parameters v and 9 satisfy the following inequalities: 

(4.23) v min < v < +oo, sC 9 sC ir. 

From Fig. 4.2 and Fig. 4.3 we see that for these values of parameters graphs of 
functions (4.21) and (4.22) intersect at unique point and determine smooth function 
z = z(v, 9) satisfying the equation (4.15). For each fixed value of v from the domain 



determined by inequalities (4.23) the function z(v, 9) is increasing in 9, it takes all 
values from closed interval [—w/2, +7r/2]: 

ITT . 7T 

0=0 1 le=7r 1 

Exactly at one point 9q = ^o(^) m th° interval [0, 7r] this function vanishes, while 
its derivative in 9 at this point is equal to unity: 

(4.25) z e (v,e (v)) = l. 

Note also that any fixed value of 9 the function z(v, 9) is increasing function in v, 
though the interval of its values here is more narrow, and it depends on 9. 

Function z(v,9), which is constructed graphically, determines the function b = 
cotanz = b(v, 9). For the fixed value of v it is decreasing function in 9, but it has 
an infinite break at the point O = 9 (v). From the equality (4.25) we derive 

(4.26) 6( U; #) = _2_ + (i) for _ Q . 

From (4.24) we obtain that b(v,9) vanishes at both ends of interval [— ir/2, +n/2]: 

(4.27) 61 =0, 61 =0. 

10=0 le=7r 
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Using b(v, 9), now we define the function A(v, 9) by the following formula 



(4.28) 



A(v, 9) = exp 



.p. J b(v. 



\ 



t) dr 



\ -*/2 



/ 



From (4.26) it follows that the function (4.28) vanishes at the point 9q = 9o(v), while 
from (4.27) we derive boundary conditions (4.19) for this function. By construction 
the function (4.28) is a solution of the equation (4.13). However its domain (4.23) 
do not embrace the whole phase space. In order to expand its domain one should 
note that if A(v,9) is the solution of the equation (4.13), then the product C(v) ■ 
A(v, 9), where C — C(v) is an arbitrary smooth function, is also the solution of this 
equation. Let's choose the function C(v) such that 

C(v) = for v ^ w min , 
C{v) = 1 for v ^ vq > t> min . 

Now, instead of formula (4.28), we define the function A(v, 9) by formula 



(4.29) 



A(v, 9) = C(v) • exp 



V 



■J b(v, 

-tt/2 



\ 



t) dr 



) 



Function (4.29) is determined in all phase space, except for those points, where 
v = |v| = 0. It is the solution of the equation (4.13) and it satisfies boundary 
conditions (4.19). 

5. Theorem on non-coincidence of classes. 

The solution of the equation (4.13) constructed by formula (4.29) determines 
the function ^(w 1 , . . . , v n ), which, in turn, determines force field F of Newtonian 
dynamical system in R™ admitting the normal blow-up of points. The construction 
of this field has functional arbitrariness due to the function / = f(w) in (4.20) and 
the function C = C(v) in (4.29). Formula (3.6) for scalar field A corresponding to 
dynamical systems admitting the normal shift of hypersurfaces also has functional 
arbitrariness due to the function W — IV^x 1 , . . . , x n , v). Let's study which part of 
this arbitrariness remains if we assume A to be spatially homogeneous function with 
axially symmetric dependence on v. The dependence on the direction of velocity 
vector v in (3.6) is completely determined by the sum 



(5.1) 



B = 



•E 

i=l 



W v 



E 



v i y,w _ (v | vw) 



w v 



w v 



Changing v by —v, we change the sign of this sum, but the value of v — |v|, which 
is the argument of function W, remains unchanged. For the case h = 1 this yields 



(5.2) 



1 

w v 



A(v)+A(-v) 
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In spatially homogeneous case right hand side of (5.2) doesn't depend on coordinates 
x 1 , ... , x n . Hence for h = 1 the quantity 1/W V depends only on v. Let's denote 
it by H(v). For the quantity B this yields 



(5.3) B = 



A(v) - H(v) for h=l, 
A(v) for h = 0. 



Thus, for both cases the quantity B = B(v) in (5.3) do not depend on coordinates 
x 1 ,. .. ,x n . Let e ij • • ■ j e n be unitary vectors directed along coordinate axes. Sub- 
stituting v = v ■ &i into the sum (5.1), we get 

(5 ' 4) m * = wT = 

The quantities rrii in (5.4) determine some vector m. Due to (5.3) they do not 
depend onx 1 , ... , x n . Hence m = m(v). Substituting the quantities (5.4) into 
the sum (5.1) and further into the formula (3.6) for scalar field A, we get 



(5.5) 



H(v) + (v | m) for h = 1, 
(v | m) for h = 0. 



Conclusion: the condition of spatial homogeneity reduces functional arbitrariness 
in (3.6) to the choice of (n + 1) functions of one variable. These are the function 
H(v) and components of the vector m(v) in formula (5.5). 

For the fixed value of v = |v| the function A(v) in (5.5) possess axial symmetry 
with the axis directed along the vector m(v). While for the function (4.29) the axis 
of symmetry doesn't depend on v, it is directed along the vector e„. Aiming to 
express the function (4.29) by formula (5.5), we should choose 



m(v) = ■ e r , 



In variables v and 9 in two-dimensional plane of axial section this yields 

, f H(v) + C(v) cos9 for h = 1, 

5.6 A = \ 

y ' 1 C(v)cos6 for/i = 0. 

Similar to (4.29), formula (5.6) contain functional arbitrariness determined by two 
functions of one variable. But this functional arbitrariness does not affect the de- 
pendence of A upon angular variable 0. While the dependence on 9 in formula 
(4.29) is much more complicated. It is determined by the choice of function f(w) 
in the equation (4.20); in general case it is not reduced to trigonometric function 
y = cos 9. Thus, formula (4.29) determines some solution of the system of weak nor- 
mality equations (1.4), which is not the solution for additional normality equations 
(1.5). Therefore we can formulate the main result of present paper. 

Theorem 5.1. In multidimensional case n ^ 3 class of Newtonian dynamical sys- 
tems admitting normal blow-up of points is the expansion of the class of systems ad- 
mitting the normal shift of hyper surf aces, and it doesn't coincide with the latter one. 
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Note that in thesis [20] for two-dimensional case n = 2 Andrey Boldin has 
constructed another (more explicit) solution for the equation (4.13). It is expressed 
through elliptic functions. However, it is local and it doesn't satisfy the conditions 
(4.19). Therefore this solution cannot be used in multidimensional case n 3 for 
proving theorem 5.1. 
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